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ABSTRACT 

■ We discuss dust properties in the interstellar medium (ISM) of nearby normal galax- 
qq ' ies, by comparing observations in the ultraviolet (UV) with simulations by a radiative 

transfer model. The observed UV colours of nearby galaxies show a reddening relative 

■ to their expected intrinsic colours. Some authors argued that the Milky Way dust 
' cannot reproduce the reddening because of the prominent 2175 A absorption bump. 
. Other authors proposed a reduction mechanism of the bump strength in an attenua- 
' tion law derived from the ratio of the observed intensity to the intrinsic one through an 

age-selective attenuation (i.e., young stars are more attenuated selectively). We newly 
find that the wavelength dependence of the scattering albedo also has a strong effect 
I 1 on the UV colour; an albedo decreasing toward shorter wavelengths (except for the 

absorption bump range) produces a significant UV reddening. After comparing the 
-j-j ■ observed UV colours of nearby normal galaxies with those expected from radiative 

transfer simulations assumed several dust models, we find two sorts of dust suitable 
for these galaxies: (1) dust with a bump and a smaller albedo for a shorter wave- 
length (except for the bump range), and (2) dust without any bump but with an 
almost constant albedo. If very small carbonaceous grains responsible for the common 
unidentified infrared emission band are also the bump carrier, the former dust is fa- 
vorable. Finally, we derive mean attenuation laws of various dust models as a function 
of the UV attenuation, and derive some relations between the UV attenuation and 
observable /theoretical quantities . 

Key words: dust, extinction — galaxies: ISM — infrared: galaxies — radiative 
transfer — ultraviolet: galaxies 
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1 INTRODUCTION 

Dust grains strongly affect our observations, especially in 
the ultra-violet (UV) band, through extinction (absorption 
and scattering processes). When we observe galaxies whose 
stars are not resolved by a telescope, the effects of the radia- 
tive transfer, i.e., multiple scatterings and the configuration 
of dust grains and stars, are also effective. In this paper, 
we call the effective extinction including such effects atten- 
uation 1 . To correctly understand the intrinsic spectrum of 
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1 Other words are also used in the literatur e: obscuration, effec- 

tive extinction, and so on. See discussions in ICalzettil feOOllV 



galaxies (unresolved into individual stars), we must correct 
observational data for the dust attenuation not for the ex- 
tinction. 

Theoretically, to know the dust attenuation through a 
galactic disc is to solve the equation of the radiative transfer 
through the disc, assuming the wavelength dependence of 
average dust properties of absorption and scattering, i.e., the 
extinction law. The most prominent feature in the average 
extinction law of the Milky Way (MW) is the absorption 
"bump" at 2175 A. We find this feature toward almost all 
lines of sight in the int erstellar medium (ISM) of the MW 
(e.g., lFitzpatricklll999l) . suggesting that the carrier of the 
bump is quite common in the ISM of the MW. We also 
find the bump in the average extinction law of the Large 
Magellanic Cloud (LMC) but not in the average extinction 
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law o f the Small Magellanic Cloud (SMC) (e.g., IWhitted 
l2003f) . 

Observationally, the wavelength dependence of the at- 
tenuation amount, i.e., the attenuation law, has been ob- 
tained for nearby UV bright starburst galaxies observed with 
the IUE satellite dCalzetti et aljll994f) . In an average atten- 
uation law of these galaxies ( so-called the Calzett i law), the 
bump is very weak or absent dCalzetti et alj Fl994*l , although 
some ga laxies in the sample s how a sign of the bump in their 
spectra dNoll fc Pierinill2005ri . Based on a radiative transfer 
model. iGordon et alJ dl99lT and lWitt fc Gordonl fcOOOt) ar- 
gued that this lack of the bump in the Calzetti law is a sign of 
the absence of the bump in the extinction law, suggesting the 
absence of the bump carrier in the starburst region. On the 
other hand, with a nother radiative transf er model, GRASIL 
dSilva et al]ll99Sl) . lGranato et all d2000l) pointed out that, 
even if the bump exists in the extinction law, the strength 
of the bump can be greatly reduced in the attenuation law 
by a radiative transfer effect coupled with an age-selective 
attenuation, i.e. young stars are more attenuated selectively. 

The IUE starburst galaxies also show a tight correlation 
between the observed UV spectral slope (/3; fx oc A"' 3 ) and 
the infrared (IR)-to-UV flux ratio (so-called IRX): a red- 
der 8 for a larger IRX dCalzetti et all 11994 iMeurer et"aH 
1999) . Since the IRX well relates to the UV attenuation 
l^uatfc XuliggH) . the correlation means that the UV spec- 
trum (i.e. the UV colour) becomes redder mon otonically 
as the attenuation increases. IWitt fc Gordonl i2000t) argued 
that the UV colour cannot be redden by the extinction law 
of the MW even if the dust column density of a medium in- 
creases. This is because the absorption bump lies in the near- 
UV (NUV). Indeed, the extinction in the NUV is slightly 
larger than that in the far-UV (FUV) for the MW extinc- 
tion law. Their finding again suggests the lack of the bump 
carrier in the starburst galaxies. 

A suggested candidate of the bump carrier is very small 
carbona ceous grains like P AHs (Polycyclic Aromatic Hydro- 
carbons; iLeeer et alll989ft QCC s (Quenched C arb onaceous 
Composites: ISakataetHul Il983l : IWada et al.l Il999t) . and 
UV pr ocessed HACs (Hydr ogenerated Amorphous Carbon 
grains -[ Mennella et al.ll998l), although this is not settled yet 
dWhittedl2003l : |Prainell2003al : iHenning. Jaeer. fc Mutschkel 
l20o3) ~These very small carbonaceous particles are confi- 
dently attributed t o the u nidentified infr ared (UIR) emission 
band in 3-13 urn dLeger fc PueelJll984l ISakata et allll984 
IWhittelil2003T) . The UIR emi ssion ban d is quite common in 
the I SM of the MW fe.g..lOnakal20oj) and of other galaxies 
(e.g., iGenzel fc Cesarskvl |200cIT except for low-metallicity 
(< 1/5 Z/Zp,) galaxies in which the UIR emission is weak or 
absent dEngelbracht et al-l koOS ) . If the very small carbona- 
ceous grains producing the UIR emission are really responsi- 
ble for the bump, we should find the bump in the extinction 
law of other galaxies (b ut not so low-metal licity). Indeed, 
the bump found in M31 dBianchi et al.lll996Tl and some dis- 
;xample, a lensing galaxy at z = 0.83 
and Mg II absorption systems at z — 1.5 
There are also signs of the bump im- 
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printed in the observed U V spectra of a galaxy at z = 0.048 
Bjn^areUa^^L|j2005a), of some IUE starburst galaxies 



Noll fc Pierinill2005D and of some star-forming galaxies at 
2 llNoll fc Pierinill2005l) . 

The amount of the bump carrier may depend on the 



star-forming activit y ||G ordorl2o'o5l . and references therein) . 
IGordon fc Clavtonl |l99^fouiid one sight line with the 
bump in the SMC. Interes tingly, this sigh t line is toward a 
quiescent area in the SMC. IValencic et alJ J2003) found one 
sight line without the bump toward an actively s tar-forming 
regio n (Trumpler 37 ) in t he MW (see also ISofia et alJ 



I2005fl . IWhittet et"aD J2004D also found a sight line with- 
out the bump towar d a molecular cloud in the MW (see 
also ISofia et al. 2005). Furthermore, the UIR emission flux 
(possibly related to the bump carrier) relative to the far-IR 
flux decreases as the intensi ty of the interstellar radiation 
field increases (IOnakall2004h . Such a bump carrier fragile 
against the star-forming activity may explain the absence of 
the bump in the Calzetti law. On the other hand, there is 
a good correlation between the UIR emission strength and 
the star-forming activity for starburst gala xies as well as for 
normal galaxies (IGenzel fc Cesarskvl l2000'). 

Quiescent or modest star-forming "normal" galaxies 
and ultra-luminous infrared galaxies (ULIRGs) do not fol- 
low the tight correl ation of the IUE starburst galax ies on 
the IRX-3 diagram jBelll 120021: fColdader et alJ 120021) . Par- 
ticularly, normal galaxies show systematically redder UV 
colours than thos e of the IUE starburst galaxies teellll2002l : 
iKong et al.ll2004h . This fa ct has recently been confirmed 
by the GALEX satellite dMartin et al.1 120051) for larger 
sam ples of nearby ga l axies selected in th e NUV or opti- 
cal dBuat et all 120051 : ISeibert et aTll2005l) . Resolved star- 
forming regions (aperture size of 520 pc) in M51 also show 
the same trend as the normal galaxies on th e IRX-UV colour 
diagra m llCalzetti et alj20 05'l. According to lWitt fc Gordonl 
(2000), the MW type dust can not reproduce even the UV 
colour of the IUE starburst galaxies, much less the redder 
UV colour of normal galaxies. Does the red UV colour of nor- 
mal galaxies indicate the lack of the bump carrier in their 
ISM and suggest different origins of the bump and the UIR 
emission? 

With a simple power-law typ e attenuation l a w (i.e . 
without a bump) as introduced bv ICharlot fc Fail j.2000), 
IKong et al.l d2004l) claimed that the redder UV colour of 
normal galaxies is due to a post-burst stellar population 
if galaxies have an intermittent star formation. They ex- 
pected a trend that the current star formation rate relative 
to the past average one (so-called the birth-rate parame- 
ter) depends on the distance from the starburst relation on 
the IRX-/3 diagram . However, the obser ved trend is weak 
dCortese et alJ 120061: Ipanuzzo et al.ll2006f) . In addition, the 
expected recent burst of nearby galaxies observed with the 
GALEX see ms to be too weak to ch ange their UV colours 
signifi cantly (iBurgarella et alJl2005bT) . 

iBurgarella et all (l2005bT) applied a more realistic atten- 
uation law, i.e., a power-law plus a Gaussian bump to their 
statistical investigation in order to understand the nature 
of GALEX galaxies. They found that an attenuation law 
with a bump and a somewhat steep slope is suitable for 
these galaxies. Such a steep a ttenuation l aw is expected from 
an age-selective attenuation dlnoudl2005l} . Mor e recently, an 
updated GRASIL model dPanuzzo et alJ 2006) very well re- 
produced the red UV colours of GALEX galaxies with the 
MW type dust. They adopted a more realistic stellar dis- 
tribution; younger stars are more deeply embedded in the 
dust disc, whereas older stars distribute m ore extensively 
(e.g., iRobin et al.ll2003l: IZaritskv et ai1l2004l) . This realistic 
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configuration of dust and stars depending on the stellar age 
produces an age-selective attenuation. This results in a steep 
attenuation law which overcomes the blueing by the bump. 

In addition to the age-selective attenuation, this pa- 
per newly discusses the effect of the wavelength dependence 
of the scattering albedo o n the UV colour. In fact, the 
dust properties adopted bv lWitt fc Gordon! l|200(t ) and the 
GRASIL are different from each other, especially the wave- 
length dependence of albedos. This point sig nificantly affects 
the ex pected UV colour as shown later fg.Tl. lWitt fc Gordonl 
( 2000) empirically derived the wavelength dependence of the 
albedo from a large compilation of the albedos e stimated 
from observations with radiative transfer models ijOordonl 
120041 . for a review). On the other hand, the dust properties 
adopted in the GRASIL are a theoretical model by Drain e 
and co-workers l)Drainel2003bl : IWeingart,ner fc Drainel200ll) . 
Such a difference in the adopted dust models could play a 
role in the contradictory conclusions from the two groups. 

This paper thoroughly examines dust properties in 
nearby galaxies, in particular the presence of the bump 
and the wavelength dependence of the albedo, based on 
the GALEX colour. We adopt a one-dime nsional plane- 
parallel radiative transfer model developed bv llnouel {2005). 
While its computational geometry is one-dimensional, this 
model can treat the dumpiness of stars and dust thanks 
to the mega-grain approximation JVarosi fc Dwekl ll99St) . 
Th e validity of this appro ximation has been clearly shown 
by IVarosi fc Dwekl lll999l) who compared the approxima- 
tion with a three-dimensional Monte Carlo radiative trans- 
fer. Owing to the computational cheapness of the one- 
dimensional calculation, we can investigate a very wide 
range of physical quantities of disc galaxies. 

In the next section, we present a description of the 
radiative transfer model, detailed explanations of the dust 
models, and the set up of the plane-parallel discs. The main 
results and discussions are presented in section 3 where we 
search suitable dust models for nearby "normal" galaxies. 
In section 4, we derive mean attenuation laws as a function 
of the UV attenuation for future use. The final section is a 
summary of our conclusions. 



2 RADIATIVE TRANSFER THROUGH A 
GALACTIC DISC 

In this paper, we use a one-dimensional plane-parallel ra- 
diative transfer model through a galactic dis c with clump y 
distributions of stars and dust developed bv llnouel {2005). 
The dumpiness of the medium (i.e. the dust distribution) 
is treated by the mega- grain approximation which was 
first proposed bv iNeufeldl (Il99llj_ and further develope d by 
iHobson fc Padmanl (11993) and IVarosi fc Dwekl dl999T) . In 
this approximation, we regard a dusty clump as a huge par- 
ticle producing absorption and scatt ering effects like a nor - 
mal single dust grain. We note that IVarosi fc Dwekl dl999T) 
clearly show the validity of the approximation by compar- 
isons between the approximate solutions and their three- 
dimensional Monte Carlo r adiative tran sfer solutions. The 
radiative transfer code of llnouel <l2005h can also treat a 
smooth medium if we do not use the mega-grain approxi- 
mation. Since the global geometry is one-dimensional plane- 



parallel, we do not consider the bulge and the radial struc- 
ture of the disc. 

By solving the radiative transfer equations in a single 
configuration of stars and dust, we obtain a transmission 
rate curve, T\ , which is the ratio of the observable and the 
intrinsic intensities as a function of the wavelength, A. When 
there are some stellar populations with different configura- 
tions relative to dust in a gala ctic disc, the total transmission 
rate through the disc is fe.g.. lTuffs et al.ll2004l) 



(1) 



where Ti,\ is the transmission rate through the i-th config- 
uration and fi t \ is the luminosity weight of the i-th stellar 
population with a normalization of - fi,\ = 1. Equation 
(1) means that we can solve the radiative transfer equa- 
tions for each stellar population (with each configuration 
relative to dust) independently. To obtain the total trans- 
mission rate, we simply sum all the transmission rates with 
luminosity weights (§2.4). 

The transmission rate of each configuration T it \ de- 
pends on 

• the dust model, 

• the ISM model, and 

• the stellar distribution. 

If we divide stars into some populations depending on their 
age, the luminosity weight f it \ is determined by 

• the age criteria, and 

• the star formation history (SFH). 

In the following, we describe the dust models consid- 
ered here (§2.1), the ISM model (§2.2), and the stellar pop- 
ulations and distributions (§2.3). Then, we describe how to 
composite these populations (§2.4). 



2.1 Dust models 

We consider six dust models in this paper. Table 1 is a 
summary of the models and their references. There are 
two origins of these m odels; one is an empirical model by 
IWitt fc Gordonl l l200Cl) ( hereafter WG dust) and the other 
is a theoretical model bv lWeingartner fc Drainel i200ll) and 
iDraind i2003bl) (hereafter Draine dust). There are also four 
types of dust compositions: the MW, the SMC, and the two 
different LMC types (LMC av and LMC 2). The LMC av 
type is an average dust composition over many sight lines 
toward the LMC, except for the supershell region around the 
30 Doradus, and the LMC 2 type is the dust composition 
toward the supershell. In the Draine dust, the bump is as- 
sumed to be produced by very small carbonaceous particles 
includ ing the "astronom ical" PAHs designed to fit observa- 
tions iLi fc Drainell200ll) . 

Fig. 1 shows the extinction cross sections (panel [a]) 
and the albedos (panel [b]) of these models as a function 
of the wavelength. The extinction cross sections (i.e. extinc- 
tion laws) are very similar between the WG dust and the 
Draine dust if we compare the same composition type. We 
find a prominent bump at 0.22 \xm (1/A = 4.6 u.m -1 ) in the 
MW and the LMC av types, a weak bump in the LMC 2 
type, and no bump in the SMC type. We note here that the 
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Table 1. Dust models 



Model 



Reference 



MW (WG) Witt fc Gordon (20001 

SMC (WG) Witt fe Gordon (2000) 
MW (D) Draine (2003b) 

SMC (D) Weingartner fe Draine (20011 

LMC av (D) Weingartner fc Draine (20011 

LMC 2 (D) Weingartner fc Draine (20011 



Visual extinction cross section per unit dust mass (cm 2 g 1 ). 



2.60 X 10 4 
1.56 x 10 4 
2.60 x 10 4 
1.56 x 10 4 
1.95 x 10 4 
1.89 x 10 4 



10 F 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 

: (a) NUV. ,FUV 



8 - 



0.0 t 



Vi iv u v ' 7 




_L 



_L 



2 4 6 

1/A [/cxitT 1 ] 



10 



Figure 1. Differences of dust models. The panel (a) shows the ex- 
tinction laws (extinction cross sections normalized by those at the 
V band) and the panel (b) shows the albedos. The solid and short- 
dashed lines are the MW and the SMC types of Witt & Gordon 

(2000) , respectively. The dotted, dot-dashed, long-dashed, and 
three-dots-dashed lines are the MW, the LMC av, the LMC 2, 
and the SMC types of Draine (2003) and Weingartner & Draine 

(2001) , respectively. The two downward arrows in the panel (a) 
show the effective wavelengths of the two GALEX filters. 



LM C av type has a rather strong bump (Fitzr 
iMisselt et aflll99Slb 

On the other hand, the albedos are very different be- 
tween the WG dust and the Draine dust (panel [b]). Ex- 
cept for the bump region, albedos of the WG dust (solid 
and dashed lines) show a flat wavelength dependence in 
2(im _1 < 1/A < 8 jtm -1 , whereas those of the Draine dust 
(other lines) show a rapid decrease toward shorter wave- 
lengths. This different wavelength dependence of the albedo 
significantly affects the UV colour as shown later (§3). 

In Fig. 2, we show comparisons between albedos of the 
MW type dusts and those estimated from observations of 
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Figure 2. Albedos estimated in the interstellar medium of the 
Milky Way. (a) The diffuse Galactic light; diamonds: Mathis 
(1973), squares: Lillie & Witt (1976), triangles: Morgan, Nandy, 
Thompson (1976), arrow (lower limit): Murthy et al. (1993), 
shaded area: Murthy fc Henry (1995), circle: Witt, Friedmann, 
& Sassen (1997), asterisk: Schiminovich et al. (2001). (b) The 
reflection nebulae; circles: Witt et al. (1982), triangles: Witt et 
al. (1992), inverse-triangle: Witt et al. (1993), squares: Gordon et 
al. (1994), diamonds: Calzetti et al. (1995), asterisks: Burgh, Mc- 
Candliss, & Feldman (2002), stars: Gibson & Nordsieck (2003). 
(c) The darck clouds; circles: Mattila (1970), triangles: Fitzger- 
ald, Stephens, & Witt (1976), squares: Laureijs, Mattile, & Schnur 
(1987), diamonds: Witt, Oliveri, & Schild (1990), asterisk: Hur- 
witz (1994), inverse-triangle: Haikala et al. (1995), stars: Lehtinen 
& Mattila (1996). The solid and dotted lines are the dust models 
for the Milky Way by Witt & Gordon (2000) and Draine (2003b), 
respectively. 



the diffuse Galactic light (panel [a]), of the reflection nebu- 
lae (panel [b] ), and of the dark clouds (panel [c]) (see also 

22). 



lGordonll2004l1 . Since the estimated values show a large dis- 
persion and both dust models are still compatible with the 
data, we can not judge which model is better. 
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Table 2. Considered physical quantities of the ISM. 



n H [cm -3 ] 



Figure 3. Phase diagram: thermal pressure— hydrogen number 
density. The dotted curves indicate the thermal equilibrium 
points in the ISM of the MW (top: Galactocentric radius of 3 
kpc, middle: 8.5 kpc, and bottom: 15 kpc) reproduced from Fig. 7 
of Wolfire et al. (2003). The two solid lines are the approximate re- 
lations of two thermally stable phases, the warm neutral medium 
(WNM) and the cold neutral medium (CNM), as expressed in 
equations (2) and (3). 



2.2 ISM models 

We consider two cases of the ISM, smooth and clumpy, in a 
plane-parallel disc. We do not consider any systematic ver- 
tical structure of the disc; the mean gas (dust) density is 
constant along the vertical axis from the equatorial plane to 
a height /id. Above this height, nothing produces absorption 
and scattering. For a clumpy medium, we have clumps dis- 
tributed randomly in the gas disc, keeping the constant gas 
density in a volume average. 

To model the clumpy mediuim we assum e a multi-ph ase 
ISM picture (e.g.. lField et alJll969|. see S2 in llnouell2005l for 
details). Assuming the thermal energy and chemical equi- 
libria in the ISM with temperatures lower than 10 4 K, we 
have two thermally stable phases (see Fig. 3) : the warm neu- 
tral medium (WNM) and the cold neutral medium (CNM). 
They are regarded as the inter-clum p med ium and clumps, 
respectively. Based on IWolfire et, alJ l)2003f) . we adopt ana- 
lytical approximations of the relations between the thermal 
pressure and the density of these two phases as 



p/k E 



nn,- 



10 4 Kcm- 3 
and 



1 cm" 



(WNM), 



(i§^r (cnm) - 



(2) 



(3) 



10 4 - 5 Kcm- 3 Vl0 3 cm 

They are shown in Fig. 3 as solid lines. By assuming a mean 
thermal pressure, we have a corresponding density contrast 
between the two phases. We also assume a mean density of 
the ISM to obtain a volume filling factor of clumps. Finally, 
the clumps are assumed to be self-gravitating in order to 
specify their radius (i.e. the Jeans length). 

In order to compare the models with observed galaxies 
and to extract the information of dust properties, we cover a 
wide range of physical quantities of disc galaxies: the mean 
ISM hydrogen density (rh = 0.5-24 cm -3 ), the half height 
of the dusty disc {ha — 50-300 pc), the dust-to-gas mass 



Quantity 



Values 



Unit 



Pth/^E 
V b 
n H c 



K cm 



10 3 , 10 3 - 5 , 10 4 
10 -3 , 10 -2 ' 5 , 10 -2 
0.5, 0.75, 1.0, 1.5, 2.0, 3.0, 
4.0, 6.0, 8.0, 12.0, 16.0, 24.0 cm" : 
50, 100, 150, 200, 250, 300 pc 



a Mean thermal pressure. 
" Dust-to-gas mass ratio. 
c Mean hydrogen density. 
^ Half height of the dusty gas disc. 



ratio (T> = 0.001-0.01), and the mean ISM thermal pressure 
(Pth/&B = 10 3 0-4 K cm -3 ), which are summarized in Ta- 
ble 2. When we keep the two phases described above, allowed 
densities are restricted between the two solid lines in Fig. 3 
for a fixed thermal pressure. Thus, we can take jih = 0.5-6.0 
cm -3 for pth/k B = 10 3 K cm -3 and n H = 2.0-24.0 cm -3 
for pth/kB = 10 4 K cm -3 . For pth/te = 10 3 5 K cm -3 , we 
can take all the values of tih listed in Table 2. Consequently, 
we have 504 sets of parameters for the clumpy ISM. 

In a smooth ISM compared with the clumpy ISM later, 
we do not need to specify the ISM thermal pressure because 
we do not make the two-phase medium. Hence, we have 216 
sets of parameters for t he smooth ISM. 

As pointed out by iDopita et"al] (120051) , the ISM pres- 
sure affects the size of H II regions and the dust temperature 
in and around these regions. Since we do not consider the 
spectral shape of the dust IR emission in this paper, we omit 
this effect. In starburst galaxies, a very high ISM pressure 
like Pth/fcB = 10 6 K cm -3 is observed (iLord et alJ Il996l : 
Heckman ct al However, the two-phase equilibrium 

is not established in such a high ISM pressure as shown in 
Fig. 3 (there is no branch of the WNM solution). Since we 
stand on the two-phase model, we restrict ourselves within 
the range of Pth/fcr j = 1Q 3 ~ 4 K cm -3 which is observed 
in the ISM of MW jMverslll97Sl) . Thus, our model is for the 
"normal" galaxies in this respect. 

The total disc optical depth at the visual band (0.55 
urn) along the normal of the disc is 



n o ( ki > v 
T V = 0.2 — — — 

V 10 4 cm 2 g- 



_\ / n H 
1 ) Ucm- 3 



h A 



150 pc A 10 



V 



■(4) 



We note that this optical depth is not observable because 
it is just proportional to the input dust column density and 
does not include the radiative transfer effect. With the vi- 
sual extinction cross section per unit dust mass, kd,v, given 
in Table 1, the considered range of the input optical depth 
along the normal axis of the discs is 0.005-25, which cor- 
responds to 0.002-5 Mq pc -2 in terms of the dust column 
density (see also §4.2.3), 

2.3 Stellar populations and their distributions 

First, we consider an age-dependent scale-height of the 
stellar distribution. In the MW, the scale-height of stars 
younger than 0.1-1 Gyr is ~ 50 pc, whereas the scale- 
height of the older stars are ~ 300 pc feinnev fc Merrifieldl 
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Table 3. Properties of the stellar populations. 



Population 



Age Layering parameter Distribution 



Young < 10 Myr 3.0 Clumpy /smooth 

Intermediate 10-300 Myr 3.0/0.5 Smooth 

Old > 300 Myr 0.5 Smooth 



Age of galaxies 
E-folding time of the SFH 



10 Gyr 
5 Gyr 



Il998l: iRobin et alJl2003l) . This is an observational fact al- 
though th e origin of this age dependent height is still con- 
trove r sial llSpitzer fc Schw arzschiidlll95l t jTinslev fc Larson! 



ll97Sl: iRanal Il99ll: iRocha-Pinto et al.l|2004D . Since the 
scale-height of the neut ral hydrogen is about 150 pc 
jBinnev fc Merrifieldll99gft , we have the layering parameters 
(ratio of the dusty disc height hd to the stellar scale-height 
h*) of ~ 3 for younger stars and of ~ 0.5 for older stars. We 
simply set these layering parameters to be constant when 
the height of the dusty disc changes in our model. 

Next, we divide the younger population into two groups. 
Observationally, the youngest stars are often associated 
with the molecular clouds. This is natural because stars 
are formed in the molecular clouds. On the other hand, 
this means very inhomogeneous distribution of the youngest 
stars. To take into account such a clumpy stellar distribu- 
tion, we bury the youngest stars into clumps in the ISM 
described in the previous subsection. Eventually, we have 
three stellar populations which are called young, intermedi- 
ate, and old stellar populations in this paper (Table 3). 

Then, we introduce the age criteria of these three stellar 
populations. The young stellar population is the population 
embedded in clumps (i.e. molecular clouds). Thus, a possible 
age criterion for this population i s a life-time of mo lecular 
clouds, for example, 10 Myr (e.g.. lBlitz fc Shulll980l) . This 
time-scale is also similar to the life-time of the latest O type 
stars. On the other hand, a comparison between the number 
of the (ultra-)compact H II regions (i.e. embedded massive 
stars) and the number of the visible O stars suggests that 
only 10-20 % of all O stars are deeply embedded in molecu- 
lar clouds JWood fc Churchwelllll989l) . To take into account 
this fact, we assume a uniform distribution of young stars in 
a clump. That is, we have young stars near the surface of the 
clump as well as young stars deeply embedded. The former 
stars would correspond to visible O stars. The intermediate 
population has a smaller scale-height than that of old stars. 
Thus, the age criterion for this population is determined by 
th e observed age-depe ndence of the scale-height. According 
to lRobin et al] (I2003T) . stars with an age of 0.1-1 Gyr have 
the smallest scale-height. Here we adopt 300 Myr for the 
criterion. This is similar to the time-scale to reach the sta- 
tionarity in the UV flux for a continuous star formation. 

Table 3 is a summary of the adopted properties of the 
stellar populations. To discuss the effect of the dumpiness 
of the young population, we consider a smooth distribution 
for this population in addition to a clumpy distribution. We 
also consider two scale heights (i.e. two layering parameters) 
for the intermediate population in order to discuss the effect 
of the small scale height of this population. 

With the age criteria described above, we have luminos- 
ity fractions of the three stellar populations if we assume a 




0.1 



1.0 



A [yU-m] 



Figure 4. Luminosity fractions of (a) the young and (b) the 
intermediate stellar populations. The young stars have an age 
less than 10 Myr, and the intermediate stars have an age between 
10 Myr and 300 Myr. The dotted line is the case of a constant 
star formation rate. The solid, dashed, and dot-dashed lines are 
the cases of an exponential star formation history with a decaying 
time-scale of 5 Gyr, 3 Gyr, and 1 Gyr, respectively. The galactic 
age of 10 Gyr is assumed for all the cases. 



SFH and a spectral energy distribution (SED) of a simple 
stellar population (SSP). Here, we assume three exponen- 
tially decaying SFHs with e-folding time-scales of 1, 3, and 
5 Gyr, and a constant SFH. Fig. 4 shows luminosity frac- 
tions for the young stars (panel [a]) and for the intermediate 
stars (panel [b]). The age of the galaxy is always assumed 
to be 10 Gyr. We have used the SEDs of the SSP with the 
Solar metallicity and a Salpeter initial mass function (0.1- 
100 Mm) based o n the Padova track used in the GRASIL 
JSilva et aljll99gl) . For an e-folding time larger than about 
3 Gyr, differences are small, especially at a wavelength less 
than 0.3 \im. Hereafter, we consider the case of the e-folding 
time-scale of 5 Gyr (solid line) as a typical case. 

We comment on a small feature at around 0.3 |xm seen 
especially in panel (b). This is probably due to the Mg I 
A2852 absorption line p roduced in old stars' atmosphere 
(e.g., iPonder et alJll998ft . This line is not prominent in the 
atmosphere of younger stars, so that the luminosity fractions 
of the young stars and the intermediate stars become rela- 
tively large. Our grid of the wavelength has only 25 points 
(see Tables C1-C6), and the grid point at A = 0.286 ixm 
is affected by the Mg I line. If our wavelength resolution 
were higher, more other features would be seen. This point 
is again discussed in §4.1. 

In computations of the radiative transfer, the emissiv- 
ity decreases exponentially along the vertical axis with a 
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Table 4. Composite models and configurations of dust and stars. 



Composite model 


CCS 


CCL 


CSS 


sss 




ISM 

V rn in o' Q^pllar Hici^n r*ii'Mr»Ti 

Scale height of intermediate stars 


Clumpy 

III 11 TY1 T~~i\r 

Small 


Clumpy 

(;1 1 1 m i""v\/ 

Large 


Clumpy 

Si m n nf~ Vi 

Small 


Smooth 

jJillUULil 

Small 




Composite model 


CCS 


CCL 


CSS 


SSS 




Young 

Intermediate 

Old 


ccl3 
csl3 
csl05 


ccl3 
csl05 
csl05 


csl3 
csl3 
csl05 


ssl3 
ssl3 
ssl05 




Configuration 


ccl3 


csl3 


csl05 


ssl3 


ssl05 


ISM 

Stellar distribution 
Layering parameter a 


Clumpy 
Clumpy 
3.0 


Clumpy 
Smooth 
3.0 


Clumpy 
Smooth 
0.5 


Smooth 
Smooth 
3.0 


Smooth 
Smooth 
0.5 



a Ratio of the height of the dusty disc to the scale height of the stellar distribution. 



scale-height h* for the smooth stellar distribution. For the 
clumpy stellar distribution, the emissivity is locally reduced 
by multiplying a factor P esc whic h is the esca pe energy frac- 
tion from a clump (see eq.[10] in llnoue) l2005l). Here, we use 
an advanced version of the model~bv lfaouer feOOot) : we take 
into account the exponentially decreasing distribution with 
a scale-height h, for stars embedded in clumps. We describe 
the formulation of this improvement in appendix A. Above 
the disc height h^, there is no stars for the clumpy case 
because of the absence of clumps. 

2.4 Composite of star/dust configurations 

Now, we composite the three stellar populations intro- 
duced in the previous subsection. This is done by equa- 
tion (1). We call the dust/star geometry for a single stel- 
lar population "configuration", and call a combination of 
three stellar populations with different configurations "com- 
posite model". Considering composite models, we have 
three points to choose the set-up: the clumpy/smooth ISM, 
the clumpy/smooth young stellar distribution, and the 
small/large scale height for the intermediate stellar popula- 
tion. Here, we consider four composite models listed in the 
top part of Table 4. First, as a standard model, we consider 
the case with the clumpy ISM, the clumpy young stellar 
distribution, and the small scale height of the intermediate 
stellar population. This case is called "CCS". To examine 
the effect of the scale height of the intermediate popula- 
tion, we consider the case with the clumpy ISM, the clumpy 
young stellar distribution, and the large scale height of the 
intermediate population. This is called "CCL" . To assess the 
effect of the young stellar distribution, we consider the case 
with the clumpy ISM, the smooth distribution of the young 
stars, and the small scale height of the intermediate stars. 
This is called "CSS" . Finally, we consider the case with the 
smooth ISM, the smooth distribution of the young stars, 
and the small scale height of the intermediate stars, which 
is useful to discuss the effect of the ISM dumpiness and is 
called "SSS". 

To produce these four composite models, we need five 
configurations of dust and stars. The bottom of Table 4 



is a summary of these configurations. The first case is the 
clumpy ISM, the clumpy stellar distribution, and the layer- 
ing parameter of 3.0 (i.e. smaller scale height of the stellar 
distribution than the height of the dusty disc) . This is called 
"ccl3" and for the young stellar population in the CCS and 
the CCL models. The next is the clumpy ISM, the smooth 
stellar distribution, and the layering parameter of 3.0. This 
is called "csl3" and for the intermediate population in the 
CCS model and the young and intermediate populations in 
the CSS model. The third is the clumpy ISM, the smooth 
stellar distribution, and the layering parameter of 0.5 (i.e. 
larger scale height of the stellar distribution than the height 
of the dusty disc). This is called "csl05" and for the inter- 
mediate and old populations in the CCL model and the old 
stellar population in the CCS and CSS models. The fourth 
is the smooth ISM, the smooth stellar distribution, and the 
layering parameter of 3.0. This is called "ssl3" and for the 
young and intermediate populations in the SSS model. The 
last one is the smooth ISM, the smooth stellar distribution, 
and the layering parameter of 0.5. This is called "ssl05" and 
for the old stellar population in the SSS model. The middle 
of Table 4 shows correspondences between the configurations 
and the stellar populations in the composite models. 2 



2 For each configuration, we consider 6 dust models (§2.1) and 
504 sets of physical quantities in a clumpy ISM or 216 sets in a 
smooth ISM (§2.2). In total, we have 3,024 cases for each csl05, 
csl3, and ccl3, and 1,296 cases for each ssl05 and ssl3. After solv- 
ing the radiative transfer equation with the number of the an- 
gular coordinate of 16, we have 48,384 transmission rate curves 
(T\, transmission rates as a function of the wavelength whose 
resolution of 25 from 0.1 ixm to 3.0 nm) for each csl05, csl3, and 
ccl3, and 20,736 for each ssl05 and ssl3. Finally, we composite 
T\ of these configurations with luminosity weights as shown in 
Fig. 4 and obtain 48,384 composite T A for each CCS, CCL, and 
CSS cases, and 20,736 for the SSS case. In total, we have 165,888 
transmission rate curves for the four composite models. 
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3 GAL EX COLOURS AND DUST PROPERTIES 

In this section, we show that the GALEX colour is strongly 
affected by dust models, especially by the presence of the 
bump and by the wavelength dependence of the albedo, 
rather than by the geometry in the radiative transfer. First, 
differences among the configurations (§3.1) and the com- 
posite models (§3.2) are shown, and then, the suitable dust 
models for the galaxies observed with the GALEX are dis- 
cussed (§3.3). We note here that the filter transmission effi- 
ciencies of the GALEX two band passes are correctly taken 
into account in all calculations. 



3.1 Comparison of configurations 

Here, we compare the UV dust attenuations among differ- 
ent configurations (the bottom part of Table 4) and dust 
models. To concentrate on effects of dust models and con- 
figurations, we consider the colour excess of the two GALEX 
bands, E(FUV - NUV) = A FUV - Anuv- Its dependence 
on the assumed SED is very weak (we assume a flat SED, 
fx cx A , for simplicity in this and the next subsection). 
Fig. 5 shows the FUV attenuation, Afuv, as a function of 
E(FUV-NUV). The panels (a) to (f) correspond to the six 
dust models described in §2.1. In each panel, the solid line is 
the locus expected from the extinction law (i.e. distant uni- 
form screen geometry) and the dashed line is that from the 
Calzetti law. When making the figure, we calculated Afuv 
and Anuv (and then E(FUV — NUV)) from the transmis- 
sion rate curves obtained from the radiative transfer calcu- 
lations with each configuration. Thus, we had 48,384 (8,064 
x 6 dusts) points each for csl05, csl3, and ccl3 (clumpy ISM 
cases) and 20,736 (3,456 x 6 dusts) points each for ssl05 and 
ssl3 (smooth ISM cases) on the E(FUV - NUV)-A FUV 
plane. Then, we divided these points into several bins in 
Afuv and obtained the maximum, minimum, and median 
values of the distribution of E(FUV — NUV) for each bin, 
each dust, and each configuration. The plotted symbols indi- 
cate the median location for five configurations: diamonds, 
squares, open circles, triangles, and filled circles are ssl05, 
ssl3, csl05, csl3, and ccl3, respectively. The bin widths are 
shown as the vertical error-bars, and the full distribution 
ranges of E(FUV ~ NUV) are shown as the horizontal error- 
bars. 

Through all the panels, we see a large variation of the 
plotted points depending on the dust models. From compar- 
isons between the MW and SMC types ([a] and [b], or [c] and 
[d]), we find the effect of the bump is very large; the most 
of the MW cases distribute within \E(FUV ~ NUV)\ < 0.3 
mag because of the bump, whereas the SMC cases can reach 
at E(FUV - NUV) > 0.5 mag. Indeed, we see "blueing" 
due to the bump in the MW (WG) case (panel [a]). The 
LMC cases ([e] and [f]) are located between the MW and 
the SMC cases. Thus, we observe the trend that the UV 
colour becomes redde r as the bump becomes weaker (e.g., 
IWitt fc Gordonll2000l) . 

We newly find that the wavelength dependence of the 
albedo also strongly affects the UV colour. From compar- 
isons between two dust models with the same extinction law 
( [a] and [c] , or [b] and [d] ) , we find that the Draine dust cases 
are systematically redder than the WG dust cases. Indeed, 
the MW (WG) case shows blueing, whereas the MW (D) 



case shows a slight reddening. While the most of the points 
of the WG dust appear in the left-hand side (bluer colour) 
of the extinction law (solid line), those of the Draine dust 
are in the right-hand side (redder colour) of the extinction 
law, except for some large Afuv cases. This difference is 
caused by the difference in albedos. As shown in Fig. 1 (b), 
the albedos of the Draine dusts decrease rapidly from optical 
to UV. On the other hand, those of the WG dusts are al- 
most flat, except for the bump domain. Relative to the WG 
dusts, the Draine dusts easily absorb the FUV photons and 
scattered out the NUV photons, so that the GALEX colour 
becomes redder, i.e. E(FUV — NUV) becomes larger. 

In each panel, we compare the different configurations. 
All configurations, except for the ccl3, show a maximum 
in Afuv and a turnover of E(FUV — NUV) in the mid- 
dle of their locus. At the largest Afuv, the minimum 
E(FUV - NUV) (or maximum for the MW [WG] case) 
reaches mag, i.e. no reddening/blueing. 3 Generally, the ob- 
served intensity along a ray is determined by the sum of two 
intensities along the ray: the intensity transmitted through 
the dusty disc and the int ensity from t he source outside the 
dusty disc (see eq. [21] in llnoud 12005). If the disc opacity 
increases, the first intensity decreases, and the observed in- 
tensity is dominated by the second intensity. In this case, 
the transmission rate (the ratio of the observed to intrinsic 
intensities) reaches an asymptotic value depending on the 
relative fraction of the second intensity in the intrinsic to- 
tal intensity, i.e., the relative amount of the source outside 
the dusty disc. The maximum Afuv is originated from the 
asymptotic value. Only the layering parameter determines 
the amount of the source outside the disc when we consider 
a single configuration. Thus, the wavelength dependence of 
the transmission rate disappears (i.e. a grey attenuation), 
and then, we have E(FUV — NUV) = at the maximum 
Afuv- The ccl3 case has no source outside the dusty disc 
because there is no clump outside the disc. Thus, this case 
does not have the maximum Afuv and shows only a weak 
turnover outside the plotted region. 

The clumpy ISM cases without the embedded stars (csl3 
and csl05) are ~ 0.1-0.3 mag bluer (redder for the MW [WG] 
case) than the smooth ISM cases (ssl3 and ssl05). This is 
because a clumpy ISM is less opaque than a smooth ISM. 
However, the case with the embedded stars (ccl3) has an 
additional local opacity due to clumps which redden (blue 
for the MW[WG] case) the UV colour, so that its colour 
excess coincides with the smooth ISM case (ssl3) before the 
turnover in E(FUV — NUV) (a small Afuv) and becomes 
redder (bluer for the MW [WG] case) than the smooth case 
after the turnover (a large Afuv)- 

In summary, the presence of the bump and the UV 
wavelength dependence of the albedo strongly affect the UV 
colour. Except for the case with the embedded stars, the 
layering parameter puts the maximum UV attenuation at 
which the colour excess becomes zero. There is also the max- 
imum (or minimum for the blueing case) colour excess. The 



3 The displacements of the symbol location from E(FUV — 
NUV) = at the largest Apuv in Fig- 5 are due to the 
binning effect. Note that the symbol location is median of the 
E(FUV - NUV) distribution in each bin. 
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Figure 5. FUV attenuation Apjjy as a function of the colour excess E(FUV — NUV). The solid lines are the relation expected from 
the extinction law (i.e. distant uniform screen). The dashed lines are the case of the Calzetti law. The open diamonds, filled squares, 
open circles, filled triangles, and filled circles are configurations of ssl05, ssl3, csl05, csl3, and ccl3, respectively (see Table 4). The vertical 
error-bars show the bin widths, the horizontal error-bars show the full widths of the distribution of E(FUV — NUV) in each bin, and 
the symbols show the median locations of the distribution. 



effects of the ISM dumpiness and the stellar dumpiness are 
not very large before the maximum colour excess. 



3.2 Comparison of composite models 

Here, we compare the GALEX colour excess, E(FUV — 
NUV), expected from the composite models summarized 
in the top part of Table 4. Fig. 6 is made by the same way 
as Fig. 5, and shows the relation between E(FUV — NUV) 



and Afuv- Each composite model is shown by each own 
symbol: SSS (squares), CSS (triangles), CCL (open circles), 
and CCS (filled circles). Six panels of (a) to (f) show six 
dust models considered here. 

As found in Fig. 5, we find that the effect of the bump 
is large from comparisons between panels ( [a] and [b] ) or ( [c] 
and [d]) and that the wavelength dependence of the albedo 
is also important from comparisons between panels ([a] and 
[c]) or ([b] and [d]), particularly for Afuv ^ 3 mag. In- 
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E{FUV-NUV) [mag] E(FUV-NUV) [mag] 



Figure 6. Same as Fig. 5, but composite cases. The filled squares, filled triangles, open circles, and filled circles are the SSS, CSS, CCL, 
and CCS models (see Table 4), respectively. 



deed, E{FUV - NUV) of the MW types ([a] and [c]) is 
about 0.5-1.0 mag bluer than that of the SMC types ([b] 
and [d]). This is the effect of the bump. On the other hand, 
E(FUV - NUV) of the Draine dusts ([c] and [d]) is about 
0.2-0.5 mag redder than that of the WG dusts ([a] and [b]). 
This is the effect of the albedo. Thus, the rapidly decreasing 
albedo toward shorter wavelengths like the Draine dust can 
partially compensate the blueing effect of the bump. 

As seen in Fig. 5, we observe a turnover in E(FUV — 
NUV) for some cases. However, the final points (the largest 
Afuv) hi each composite model do not reach E(FUV — 
NUV) = 0. This shift from E{FUV - NUV) = is caused 



by the composite process, i.e. the effect of the age-dependent 
attenuation. This effect can be understood by equation (1). 
If the disc opacity is large enough, the attenuations for any 
populations become grey; Ti t \ loses the wavelength depen- 
dence. Even in this case, we still have the wavelength depen- 
dence of fi t \. Thus, E(FUV — NUV) for an opaque disc is 
determined by the wavelength dependence of the luminosity 
weight, in other words, the SFH and the age criteria of the 
stellar populations. More interestingly, the attenuation law 
becomes independent of the dust model in such a case. In- 
deed, we find that the locations of the largest Afuv point 
of each composite model in six panels (i.e. six dust mod- 
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Figure 7. Same as Fig. 5, but the CCS model (Table 4) with different star formation histories. The open circles, filled triangles, and 
filled circles are the cases with e-folding time-scales of 1 Gyr, 3 Gyr, and 5 Gyr, respectively. The filled squares are the cases with a 
constant star formation rate. 



els) are very similar to each other. This causes rather large 
E(FUV - NUV) even for the MW dust (panels [a] and [c]). 
We can expect a red GALE X colour even with the bump 
(see also IPanuzzo et, afl l^OOfi') . The age-dependent attenua- 
tion reduces the effect of the bump through the wavelength 
dependence of the luminosity weights. 

In each panel, when we compare the CCS and CCL 
models, the effect of the intermediate population can be un- 
derstood. With a large scale-height of the intermediate pop- 
ulation (CCL), we can not reach Afuv > 2.5 mag, whereas 
we reach Afuv ~ 5 mag with a small scale-height of the 



population. As shown by iBuat et alJ {2005), FUV attenu- 
ations of many nearby galaxies selected by FIR are more 
than 2.5 mag. Thus, an intermediate population with a small 
scale-height is likely to be required. When we compare the 
CCS and CSS or SSS models, we find the effect of the clumpy 
young stars. Because of the additional local opacity due to 
clumps, the CCS model reaches ~ 1 mag larger Afuv than 
the CSS and SSS models. We find that the effect of the ISM 
dumpiness is small for Afuv % 3 mag from a comparison 
between the CSS and SSS models. In the following discus- 
sions, we deal with only the CCS model which consists of 
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the young population embedded in clumps, the intermediate 
population with a small-scale height, and the old population 
with a large-scale height, i.e. the most realistic model in this 
paper. 

Fig. 7 shows the effect of the SFH on the E(FUV - 
NUV)-Afuv relation for the CCS model. For an e-folding 
time-scale longer than 3 Gyr, the difference in E(FUV — 
NUV) is almost zero for Afuv < 2 mag, and small shifts 
(~ 0.1 mag) of the largest Afuv are observed. This small 
difference is due to the small difference of the luminosity 
weights in the UV (A < 0.3 \im) as shown in Fig. 4. On the 
other hand, a short e-folding time-scale like 1 Gyr produces 
a difference on the E(FUV — NUV)-Afuv plane. This case 
reaches only Afuv < 2 mag because the radiation is dom- 
inated by the old population emitting from the outside of 
the dusty disc. 

In summary, the effect on the GALEX colour signifi- 
cantly depends on the dust model (with/without the bump 
and the albedo) for a small UV attenuation, whereas it be- 
comes independent of the dust model for the largest UV 
attenuation. The clumpy young stars and the intermediate 
population with a small scale-height increase the UV atten- 
uation. The effect of the intermediate population is larger 
than that of the young population, at least in the SFH con- 
sidered here. On the other hand, the ISM dumpiness has a 
relatively small effect on the UV colour excess. The SFH has 
a very small effect on the UV colour excess if we assume a 
smooth SFH with an e-folding time-scale larger than 3 Gyr. 



3.3 IR-to-UV flux ratio and GALEX colour 

As shown in the previous subsection, the GALEX colour 
is very sensitive to the presence of the bump and the 
wavelength dependence of the albedo. Conversely, we may 
assess the dust models by comparing with the observed 
GALEX colours. Fig. 8 shows the diagram of the IRX 
(dust IR-to-UV flux ratio, F dnst /F FUV ) and the GALEX 
colour, FUV — NUV. The crosses and diamonds are the 
observed data of the nearby g alaxies selec t ed by NUV and 
FIR, respectively, taken fr om iBuat et alJ (120051) (see also 
llglesias-Paramo et al.l2006T) . The filled circles and error-bars 
are made by the same procedure done in Fig. 5. We show 
only the CCS model with the e-folding time-scale of 5 Gyr 
for the comparison. We calculated the model dust IR flux 
as the total absorbed flux. This procedure is described in 
appendix B. The UV flux is defined as the UV flux density 
multiplied by the effective wavelength of the UV filter. The 
solid and dashed lines in Fig. 8 are the loci expected from 
the extinction law and the Calzetti law, respectively. 

By comparing the model points with the observed data, 
we assess which dust model is suitable for the observed 
galaxies. The SMC (WG) case shows a very good agree- 
ment with the data of the UV selected galaxies and the FIR 
selected galaxies with Jd us t /Ffuv i$ 100 (panel [b]). The 
LMC av (D) and the LMC 2 (D) cases are also compatible 
with the data (panels [e] and [f]). The colours predicted by 
the MW (D) case are still ~ 0.2-0.3 mag bluer than the 
observed ones because of a strong bump and a shallow UV 
slope in the extinction law (panel [c]), although a differ- 
ent IMF like a Kroupa IMF could reduce the discrepancy 



iPanuzzo et aljfeoOrf) . 4 On the other hand, the predicted 
colours of the MW (WG) case are largely separated from the 
observed data, say ~ 0.5 mag at F& us t/ Ffuv ~ 10 (panel 
[a]). For the SMC (D) case (panel [d]), the predicted colours 
are too red (~ 0.5 mag at Faust I Ffuv ~ 10) because of a 
rapid decrease of the albedo between the two GALEX bands 
as shown in Fig. 1 (b). 

Interestingly, there is no model which reproduces the 
FIR selected galaxies w ith Jduat/ Ffuv > 100 where 
ULIRGs also distribute dGoldader et al.1 12002). 5 As dis- 
cussed in §3.2, for an opaque disc, we expect to have an 
attenuation law (i.e. transmission rate curve) independent 
of dust properties. In fact, we find that the locations of the 
most opaque point in each panel are very similar; all dust 
models predict a very similar position on the diagram for 
Fdust/ Ffuv > 100. However, the real galaxies show a very 
larg e dispersion in th e region . 

iBurgarella et all J2005b) could not explain the same 
galaxies by their analysis either and suggested that an effect 
of "decoupling" is important; the stellar populations produc- 
ing the UV and the IR are completely different. For exam- 
ple, the UV radiation comes from the population outside the 
obscured region, whereas the population heating dust which 
emits the IR radiation is embedded there. In this case, the 
UV colour is decoupled with the UV attenuation traced by 
Fdust /Ffuv- 6 In the framework of this paper, such a "de- 
coupling" would take place if we consider an intermittent 
SFH. Now, we have three stellar populations: young and in- 
termediate ones embedded in clumps and in the dusty disc, 
and old one distributed diffusely to the outside of the disc. 
Under an intermittent SFH with a time-scale longer than 
~ 300 Myr (age threshold between the intermediate and 
old populations), we can expect that the luminosity weights 
strongly vary along the time, and then, the position of the 
most opaque case on the IRX-UV colour diagram would 
vary. 

In summary, the GALEX data are consistent with WG's 
SMC dust and Draine's dusts with the bump (LMC av, 
LMC 2, and possibly MW). If we consider that the car- 
rier of t he UIR very common am ong the nearby galax- 
ies (e.g., iGenzel fc Cesarskvl |2000T) is also the bump car- 
rier as suggested by some laboratory investigations (e.g., 
ISakata et allll983fi . dust with the bump is more favorable. 
For a flat wavelength dependence of the albedo like the WG 
dust, a strong bump seen in the MW makes the GALEX 
colour too blue. However, a weak bump like the LMC 2 



4 iPanuzzo et alJ J2006t) have reproduced the data of the NUV se- 
lected galaxies with a Draine's MW dust (but a different version 
from this paper) better than here. This may be due to some dif- 
ferences between the two papers; they treat a 2-D disc + bulge 
with a smooth medium, whereas we treat a 1-D plane- parallel disc 
with a clumpy medium. Moreover, a differ ent age criterion of the 
intermediate population may play a role; IPanuzzo et al.l J2006fl 
show that a shorter criterion gives a redder colour. 

5 Onl y a few galaxies in the FIR selected sample of lBuat et alJ 
l200Eft are ULIRGs. 

6 We can still rely on the IR-to-UV flux ratio for an indicator 
of the total UV attenuation because the flux ratio is the ratio 
of the absorbed to observed radiation energies (escaped from the 
obscured region + emitted from the outside population) based on 
the energy conservation. 
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Figure 8. Dust IR-to-FUV flux ratio and the GALEX colour. The crosses and dia monds are observed data of the NUV selected and 
the FIR selected nearby galaxies, respectively, taken from Buat et al. (2005) (see also llgle sias-Para mo et al.l200fifl . The solid and dashed 
lines correspond to the extinction law and the Calzetti law, respectively. We show the CCS model (Table 4). The points predicted from 
the model are divided into several bins in the flux ratio. The vertical error-bars are the bin width, the horizontal error-bars are the full 
width of the distribution of FUV — NUV in each bin, and the symbols are the median of the distribution. 



type could be still compatible with the GALEX data. Some 
FIR selected galaxies which have a large IR-to-UV flux ratio 
(-Fdust /Ffuv > 100) cannot be explained by a smooth SFH 
model, may suggesting their recent episodic star-formation. 

4 MEAN ATTENUATION LAWS 

From the calculated transmission rate curves, we derive 
mean attenuation laws which would be useful for correcting 



the observed data of galaxies for the dust attenuation and 
for predicting the observable SEDs in theoretical works. We 
show only the CCS model (Table 4) in this section. 



4.1 Mean UV-to-NIR attenuation laws 

We find that the FUV attenuation, Afuv, represents the 
global shape of the attenuation laws. In Fig. 9, the colour 
excesses relative to Afuv at various wavelengths, E(X — 
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Figure 9. Colour excesses relative to the FUV at various wavelengths, E(X — FUV), as a function of the FUV attenuation, Apuv ■ 
The filled circles, filled triangles, filled squares, open circles, open triangles, and open squares are the colour excesses at A = 0.100, 
0.157, 0.229, 0.300, 0.599, and 1.19 um, respectively. The horizontal error-bars indicate the Apjjy bin width. The plotted points and 
the vertical error-bars are the mean and the standard deviation of the E(\ — FUV) distribution in each bin. The thin solid curves are 
the 3rd order polynomial fits. 



FUV) = A\ — Afuv, are shown as a function of Afuv- To 
make the figure, we first averaged the transmission rates, T\, 
over the angle between a ray and the disc normal (i.e. incli- 
nation angle) via equation (B2). Then, we divided the incli- 
nation averaged T\ into some bins in Afuv and calculated 
the mean and the standard deviation in the distribution of 
E(X — FUV) in each bin and each dust model. The location 
of symbols and the vertical error-bars in Fig. 9 show the 
mean (not the median) and the standard deviation (not the 



full width), respectively. The horizontal error-bars indicate 
the bin widths. Fig. 9 shows that E(X — FUV) can be ex- 
pressed as a function of Afuv very nicely. Indeed, a typical 
standard deviation is as small as 0.1 mag. Thus, we tried to 
fit the calculated E(X — FUV) by a 3rd order polynomial 
function of Afuv for each dust model as 



E(X - FUV) = a(X)A F uv + P(X)A F uv 2 + jWAfuv 3 .(5) 



Ultraviolet dust attenuation in galaxies 15 



CO o 



< 2 



1 - 



&0 o 

cd o 



. 2 



. 2 




4 - 



1 - 



4 - 



1 - 



- 



I I I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I 



_L 



_L 



_L 



I I I I I I I I I I I I I I I I I I I 




I I I I I I I I I I I I I I I I I I I 




_L 



_L 



2 4 6 
1/A [/im -1 ] 



2 4 6 
1/A [yum -1 ] 



10 



Figure 10. Mean attenuation laws. In each panel, eight attenuation laws are shown for Apjjy = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 
mag from the bottom to the top. 



Note that we have assumed E(\ — FUV) — when Apuv = 
0. The order of the polynomial function was determined, 
based on the Akaike's Information Criterion jTakeuchil 2000; 
iTakeuchi. Yoshikawa. fc Ishil 2000'l ; we tried to fit up to 6th 
order and confirmed that 3rd order is enough. The resulting 
parameters for each wavelength and each dust model are 
tabulated in appendix C. In Fig. 9, the fitting results are 
shown as thin solid curves. 

Fig. 10 shows the mean attenuation laws obtained by 
the polynomial fit of E(X — FUV). We observe that the 
significance of the bump reduces as Afuv increases (panels 
[a], [c], [e], and [f]). We also find a different slope in the UV 



(1/A ~ 5-7 u.m _1 ) between the two SMC cases (panels [b] 
and [d]); a steeper rise of the SMC (D) case is due to a rapid 
decline of its albedo (Fig. 1 [b]). We also note that the shape 
of the attenuation law for the most opaque case (top curves, 
Afuv = 4 mag) in each panel is very similar to those in 
other panels as discussed in §3.2. 

We comment on the small feature seen at 1/A ~ 3.5 
|t.m _1 in all panels of Fig. 10. This is due to the small fea- 
ture seen in Fig. 4, the stellar Mg I A2852 absorption line. 
This line is prominent in old stars' atmosphere but not in 
young and intermediate stars' one. Such an age-dependent 
strength of the line appears as a small feature in the luminos- 
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ity fractions of young and intermediate stars. As discussed 
in §3.2, the wavelength dependence of the attenuation law 
(i.e. transmission rate T\) is determined by the luminosity 
fractions if the disc is enough opaque. Thus, such a stellar 
feature appears in the attenuation law of an opaque disc. In- 
deed, the feature in Fig. 10 is clearer for larger Afuv cases. 
Interestingly, an age-dependent stellar feature can appear in 
the attenuation law through an age-selective attenuation. 

To quantify the bump reduction and the slope change 
along the FUV attenuation seen in Fig. 10, we tried to fit 
the normalized mean attenuation laws by a power law plus 
a Gaus sian bump with the same way as iBurgarella et all 
J2005bl) : 



+ B exp ( - 



A — A t 



(6) 



where the bump position Ab = 0.2175 \im and the bump 
width a = 0.02 urn. Although this function does not give 
a very good fit, we adopt it for simplicity. We left two free 
parameters: the power law index p and the bump ampli- 
tude B. Fig. 11 shows results of the fit. We find clearly that 
the global slope of the attenuation law reduces as Afuv 
increases (panel [al). This is consistent with previous in- 
vestigations (e.g . ; |VarosL^r>we3ll999|; |Feiir^ja^_et^_ajJLl999|; 
[ Witt fc Gordon! l2000t iPierini et al] 12004 iPanuzzo et alJ 
2006) For each dust model, the power law index decreases 
from a larger value than its extinction law (i.e. steeper at- 
tenuation) to a smaller value (i.e. greyer attenuation), and 
finally, the indices converge a value similar to th at of the 
Calzetti law (p = 0.7 see also lCharlot fc FaljEoOol) . 

We also find a clear trend that the bump amplitude re- 
duces as Afuv increases (panel [b]); from a similar value to 
that in the extinction law to a very small amplitude even 
less than that of the LMC 2 (D) case. Loci of the bump am- 
plitudes along Afuv are very similar, except for the LMC 2 
(D) case whose amplitude is very small even in its extinction 
law. This trend is also consistent w i th th e literature (e.g., 
I Varqsi fc Dwekl ll99St iFerrara et alJ ll99St IWitt fc Gordon! 



hood IPierini et alJl2004 IPanuzzo et al.l l200d) . The reduc- 
tion of the bump amplitude is originated from the nature 
of t he grey attenuation in an opaque disc as discussed in 
S3. 1. IBurgarella et al.l j2005bl) did not find the reduction of 
the bump amplitude along the FUV attenuation from their 
analysis of the galaxies observed with the GALEX. This 
may be because their sample consists of two different pop- 
ulations. Indeed, the bump amplitudes estimated in their 
analysis show a bimodal distribution. 

From Fig. 11, we understand why some dust models can 
reproduce the GALEX data in Fig. 8 and other can not. The 
slope of the attenuation law of the SMC (WG) is very good 
for the data. Suppose this as the reference case. Although 
the LMC av (D) and the LMC 2 (D) give a steeper slope than 
that of the SMC (WG), the bump which shifts the GALEX 
colour somewhat bluewards compensates the steepness. The 
slopes of the two MW dust cases are shallower than the 
reference SMC (WG) case, and moreover, there is the bump. 
For the SMC (D) which does not have the bump, the slope 
of the attenuation law is too steep. 
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Figure 11. Properties of mean attenuation laws: (a) global slope 
of the attenuation law and (b) amplitude of the bump. The ver- 
tical axes are defined in equation (6): index = — p and bump 
amplitude = B. The crosses and circles are MW (WG) and SMC 
(WG), respectively. The triangles, squares, diamonds, and aster- 
isks are MW (D), LMC av (D), LMC 2 (D), and SMC (D), re- 
spectively. For a reference, the parameters of the extinction laws 
and the Calzetti law are marked in the right-hand side in each 
panel. Two marks for each extinction law mean two MW laws of 
WG and Drain dusts, two LMC laws of the Draine dust, and two 
SMC laws of WG and Draine dusts. 



4.2 Relations between UV attenuation and other 
quantities 

We have derived the mean attenuation laws as a function of 
the FUV attenuation, Afuv- If we have Afuv, thus, we can 
obtain an attenuation law from the UV to the NIR. Now, we 
need to know relations between Afuv and other observable 
or theoretical quantities in order to estimate Afuv- 



4-2.1 IR-to-UV flux ratio 

The IR-to-UV flux ratio is a very go od measure of Afuv 
jBuat fc Xulll996l : iGordon et alJl2000l) . This is because the 
relation is based on the energy conservation; the IR flux 
is the flux absorbed by dust grains. Thus, the relation is 
very robust against differences of the configuration of stars 
and dust and the dust model. Indeed, Fig. 12 (a) shows 
the robustness of the relation against differences among 
dust models. The standard deviations are less than 0.2 
mag. Since the FUV attenuation is defined as Afuv = 
2.5\og{FF uv /FFUv) with F'^uv and F F uv being the in- 
trinsic and observed FUV fluxes , respectively, we c an ex- 
press the relation as follows (e.g- lMeurer et al.lll999t) : 
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Figure 12. Mean relations between the IR-to-UV flux ratio and 
(a) the FUV attenuation and (b) the fraction of the FUV ab- 
sorbed energy in the IR flux. The symbols, which almost overlap 
in the panel (a), indicate the mean relations for each dust model; 
crosses, circles, triangles, asterisks, squares, and diamonds are 
MW (WG), SMC (WG), MW (D), SMC (D), LMC av (D), and 
LMC 2 (D), respectively. The solid lines are the best fit curves 
expressed in equations (7) and (9). The dashed line in the panel 
(a) is the calibration proposed by Buat et al. (2005). 
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Figure 13. Mean relation between the FUV attenuation and the 
UV colour: (a) GALEX colour and (b) UV spectral slope 0. The 
symbols are the same as Fig. 12. 



(2005) (dashed curve in Fig. 12 [a]) is very small, except 
for the region of a small -Fdust/ Ffuv ■ Importantly, we can 
obtain a UV-to-NIR attenuation law from the observed IR- 
to-UV flux ratio with the calibration obtained here. 



A F uv = 2.h\og(l + i F vv^-\ , (7) 

V CFUVJ 

where 

e F'fuv - F FU y 

K,fuv = p , (a) 

-Fdust 

which is the energy fraction of the absorbed FUV flux in 
the IR flux. From Fig. 12 (b), we find that £, F uv is 0.3-0.6 
although the value depends on the dust model and the dis- 
persion is somewhat large. In an extremely dust poor case 
(-Fdust/-FW -> 0), Fpyy = F FUV but still F dust > be- 
cause dust grains can be exposed by ionizing photons. Thus, 
£,fuv — * when Fd ust / Ffuv — * 0. In the opposite limit, all 
stellar radiation is absorbed by dust and re-emitted in the 
IR (Ffuv -* and F dust — > F tota i, where F tota i is the total 
flux); £fuv -»• i'Fuv = F'puv / Ftotai- If we assume the SED 
of a smooth exponential SFH with an e-folding time of 5 
Gyr and a galactic age of 10 Gyr, we have ^fuv = 0.245. 
With these limits, we have found a fitting function as 

too X 

(.fuv = UV + 0.400 exp(-(0.591a; - 0.185) 2 ) , (9) 
1 + e x 

where x — log(Fd us t/fWy). This is shown as the solid curve 
in Fig. 12 (b). The relation obtained from equations (7) and 
(9) is also shown as solid curve in Fig. 12 (a). The d iffcr- 
ence of the calibration obtained here from that of lBuat et alJ 



4.2.2 UV colour 

As found bv lMeurer et al] dl999l) . the UV bright starburst 
galaxies show a tight correlation between the UV spectral 
slope (i.e. UV colour) and the IR-to-UV flux ratio (i.e. 
Afuv). However, normal galaxies do not follow the relation 
jBellll2002t Ikong et aLlliool iBuat et"al]|2005l) . Figs. 5 and 
6 show that the relation between Afuv and the GALEX 
colour excess strongly depends on the dust model. Thus, 
we expect that the UV colour and slope also depend on 
the dust model significantly. We estimated the UV slope (3, 
where fxocX 13 , from the expected flux densities at 10 wave- 
length points betwee n 0.12 and 0.26 u.m listed in table 2 of 
ICalzetti et"ai]||l994l) . As shown in Fig. 13, we find that rela- 
tions between Afuv and UV colour (or slope) show a large 
dispersion depending on the dust model. For example, Afuv 
distributes from 0.5 to 4.0 mag for FUV - NUV = 0.5 mag 
or /3 = — 1. Although we may expect a smaller dispersion of 
Afuv for a fixed UV colour (or slope) in each dust model, 
the dispersion is still ~ 1 mag. 

4-2.3 Dust column density 

To predict the observable fluxes of galaxies by theoretical 
models, we need a relation between Afuv and a theoreti- 
cal quantity. Since the total dust column density, Edust, of 
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Figure 14. Mean relation between the FUV attenuation and the 
dust column density. The symbols are the same as Fig. 12. The 
solid line is the best fit curve described in equation (10). 

Table 5. Parameters for equation (10). 
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the disc can be calculated by a galactic ch emical and dust 
amount evolution model (e.g., Ilnouel 120031) . a relation be- 
tween Afuv and Edust would be useful. Fig. 14 shows the 
^4.F!7V - £dust relation which is very robust against differences 
among dust models like the relation between Afuv and the 
IR-to-UV flux ratio. A typical standard deviation is as small 
as 0.2 mag. As Ed us t increases, first Afuv increases linearly, 
and then, becomes saturated at a certain value determined 
by the amount of the source outside of the disc, i.e. layering 
parameter (see §3.1). Taking into account such a behavior, 
we fit the data by the following function: 

3 

Afuv = ^ cn{l - exp(-biE du st)} . (10) 

i=l 

The resulting parameters are tabulated in Table 5, and the 
best fit curve is shown as the solid line in Fig. 14. 



5 CONCLUSION 

We discussed dust properties in the ISM of nearby normal 
galaxies observed through the GALEX FUV and NUV fil- 
ters. To extract the dust properties from the GALEX data, 
we built a huge set of UV-to-NIR transmission rate curves of 
disc galaxies. The set covers a very wide range of the physical 
parameters of these galaxies. First, we examined various ef- 
fects on the observed UV colour, such as the dust models, the 
age-selective attenuation including the dumpiness of stars, 
the dumpiness of the dust distribution, and the SFH. Next, 
we compared the distribution of the observational data of 
nearby normal galaxies on the IR-to-UV flux ratio and the 
GALEX colour diagram with that expected from our trans- 



mission rate curves. Then, we derived UV-to-NIR mean at- 
tenuation laws as a function of the FUV attenuation from 
our transmission rate curves. We also derived relations be- 
tween the FUV attenuation and the IR-to-UV flux ratio, the 
UV colour, and the dust column density of the disc in order 
to use the mean attenuation laws practically. 

From these analyses, we confirmed the following 
conclusions seen in the literature: the 2175 A absorp- 
tion bump, which is located in the GAL EX NUV fil- 
ter, makes UV colours m uch bluer (e.g., iGordon et alJ 
Il997t IWitt fc Gordon! l2000l) . the age-selective dust atten- 
uation (i.e. younger stars are more atten uated) reduces 
the bump strength in the attenuation law iGranato et alJ 
l200Ct iPanuzzo et alJ l200rJl , the shape of the attenuation 
law changes from steep to shallow as the FUV attenua- 
tion increases fe.g.. lWitt fc Gordonl2000l) . and the IR-to-UV 
flux ratio is a very good indicator of the FUV attenuation, 
whereas the UV colour is not good one (e.g., iBuat fc XiJ 
ll996HGordon et alilSoOch . 

Our new findings are 

(1) a rapid decline of the albedo toward a short wave- 
length (except for the bump range) makes UV colours sig- 
nificantly redder (§3.1, 3.2), 

(2) a smooth exponential SFH with an e-folding time- 
scale larger than 3 Gyr does not affect the UV colour excess 
significantly (§3.2), 

(3) the FUV attenuation is a very good measure of the 
attenuation at an arbitrary wavelength from the UV to the 
NIR in our models (§4.1), 

(4) an age-dependent stellar feature like Mg I A2752 can 
appear in the attenuation law through the age-selective at- 
tenuation (§4.1), and 

(5) the dust column density shows a very tight relation 
with the FUV attenuation (§4.3). 

From a comparison between our model and the GALEX 
observations in §3.3, we found that Witt & Gordon's SMC 
type dust shows a very good agreement with the GALEX 
data, whereas Witt & Gordon's MW type dust is not con- 
sistent with the data, and that Draine's LMC average and 
LMC 2 type dusts are also well compatible with the GALEX 
data, whereas Draine's SMC type dust is not consistent with 
the data. Draine's MW type dust gives slightly bluer colour 
than the data although this discrepancy could be resolved 
ea sily, for example, with another proper IMF as suggested 
bv lPanuzzo et all i2006t) . Therefore, the main conclusion of 
this paper is that 

(6) in the ISM of the nearby normal galaxies, there is 
either dust with a bump and a smaller albedo for a shorter 
wavelength (except for the bump range), or dust without 
any bump but with an almost constant albedo. 

If we regard very small carbonaceous grains respo nsible for 
the U IR emission band as the bump carrier (e.g., IWhittetl 
2003), the former dust is more suitable for nearby normal 
galaxies, beca use the UIR emission is quite ubiquitous in 
these galaxies <|Oenzel fc Cesarskvll2000l) . 

Explaining the behavior of the galaxies with a large IR- 
to-UV flux ratio (typically > 100) on the IR-to-UV flux ratio 
and the GALEX colour diagram remains a future work. Our 
model predicts the convergence to a certain GALEX colour 
independent of the dust model for galaxies with such a large 
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IR-to-UV flux ratio, whereas the nature shows a very large 
dispersion of the colour. This may suggest a recent episodic 
star formation, which decouples the stellar populations ob- 
served in the UV and in the IR, in these galaxies. 
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APPENDIX A: CLUMPY STELLAR 
EMISSIVITY 

Here we describe how to manage dumpiness of the stellar 
emissivity. We do not consider any systematic distribution 
of clumps, i.e. clumps distribute uniformly in the gas+dust 
disc. However, we consider that a mean number density of 
embedded stars decreases along the vertical axis of the disc. 
We update the trea tment of the clumpy stellar emissivity of 
the section 3.4.1 in llnoud 1120051) as follows. 

First, we normalize the intrinsic emissivity, r/,, as unity. 
When clumps exist only within the gas+dust disc of the half 
height hd and the embedded stars exist only in clumps, this 
normalization becomes 

rj*(z)dz = l. (Al) 

If the number distribution of the embedded stars is an ex- 
ponential function along the vertical axis, z, with an e- 
folding scale of h,, we have ?7*(.z) = exp(— \z\/h*)/2h»(l — 
exp(— h&/h*)) for \z\ < /i d , otherwise r/ t — 0. 

A part of the radiation from the embedded stars is lo- 
cally absorbed by dust in the clump where the stars are em- 
bedded. If the embedded stars distribute uniformly in the 
clump and the scattering by dust is isotropic, the photon 
escape probability from the clump, P csc (r c i, oj c i), with the 
optical depth radi us t c i and the effective albedo u; c i , is given 
by equation (10) in llnoue| i2005l) . Therefore, the clumpy stel- 
lar emissivity input into the transfer equation is 

f P °sc(t- c1 ,^ c1 ) cxp(-|z|/h») ( r I I < , \ 

n?(z) = \ 2h, { i- cxp( -h d /h,)} ( tor M «S M _ (A2 ) 

I (for \z\ > h d ) 

In this case, we do not have any incident radiation at the 
top of the disc as the upper boundary condition of the cal- 
culation. In addition, we set a mirror boundary condition at 
the equatorial plane of the disc. 



APPENDIX B: ESTIMATION OF THE DUST 
INFRARED LUMINOSITY 

Here, we explain the estimation procedure of the dust 
IR luminosity, in other words, the total absorbed lu- 
minosity. The dust IR luminosit y is expressed as (e.g., 
Ilnoue. Hirashita. fc KamavallioOoTl 

L dust = L LyQ + (1 - /lc)£lc + / La(1 - (T\))d\ , (Bl) 

where I/Lya and Llc are the Ly a emission line and the Ly- 
man continuum luminosities, L\ is the intrinsic luminosity 
density of all stellar populations, /lc is the luminosity frac- 
tion of the Lyman continuum absorbed by hydrogen atoms, 
and (T\) is the transmission rate averaged over the angle 
between the disc normal and a ray. That is, 

(Tx) = f T\(p)d[i , (B2) 
Jo 

H is the cosine of the angle. Note that we should inte- 
grate the composite transmission rate (with an angle de- 
pendence) over the angle to obtain the total absorbed en- 
ergy within the disc. We assume the same SED as that 
used to obtain the luminosity weights in equation (1) and 
Fig. 4 (solid lines). The first two terms in equation (Bl) 
can be reduced to (1 — q/lcO-^ lc with a is a numerical fac- 
tor o f the order of unity (e.g.. Ilnoue. Hirashita. fc Kamaval 
2000). With the case B approximation and the assumed 
SED ; we have a ~ 0.6. Some estimates suggest /lc — 
0.5 (Inouc. Hirashi ta. fc Kamaval 1200 ll : llnoud boOll 120021: 
iHirashita. Buat. fc Inouei2003ri . Thus, we adopt 1 — o/lc = 
0.7. This choice does not affect the estimated dust IR lumi- 
nosity significantly because a typical luminosity fraction of 
Llc in the total luminosity is small enough (15%). 



APPENDIX C: PARAMETERS FOR MEAN 
ATTENUATION LAWS 

We summarize the fitting parameters for equation (5) of each 
dust model in Tables CI to C6. 
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Table CI. Parameters for MW (WG) Table C3. Parameters for MW (D) 
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Table C2. Parameters for SMC (WG) Table C4. Parameters for LMC av (D) 
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1.333E-02 


2.383 


-9.178E-01 


-3.951E-02 


1.089E-02 


2.383 


-9.094E-01 


-4.650E-02 


1.284E-02 


3.000 


-9.379E-01 


-3.338E-02 


8.090E-03 


3.000 


-9.141E-01 


-4.493E-02 


1.214E-02 



Table C5. Parameters for LMC 2 (D) 
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7 



0.100 
0.114 
0.129 
0.143 
0.157 
0.171 
0.186 
0.200 
0.214 
0.229 
0.243 
0.257 
0.271 
0.286 
0.300 
0.378 
0.475 
0.599 
0.754 
0.949 
1.194 
1.504 
1.893 
2.383 
3.000 



9.849E-01 
5.809E-01 
3.089E-01 
1.549E-01 
-4.878E-02 
-2.238E-01 
-3.441E-01 
-3.839E-01 
-3.658E-01 
-4.411E-01 
-5.344E-01 
-5.998E-01 
-6.490E-01 
-6.673E-01 
-7.130E-01 
-8.037E-01 
-8.581E-01 
-8.795E-01 
-8.883E-01 
-8.928E-01 
-8.973E-01 
-9.019E-01 
-9.041E-01 
-9.065E-01 
-9.120E-01 



-2.498E-01 
-1.278E-01 
-5.578E-02 
-2.023E-02 
1.837E-02 
4.130E-02 
5.403E-02 
6.664E-02 
7.510E-02 
7.418E-02 
5.563E-02 
4.137E-02 
4.151E-02 
2.952E-02 
4.771E-02 
2.841E-02 
7.700E-03 
-1.402E-02 
-2.872E-02 
-3.802E-02 
-4.311E-02 
-4.537E-02 
-4.566E-02 
-4.516E-02 
-4.371E-02 



1.699E-02 
6.279E-03 
3.272E-04 
-1.482E-03 
-2.057E-03 
-5.816E-04 
4.379E-04 
-2.391E-03 
-7.153E-03 
-4.840E-03 
3.301E-03 
7.775E-03 
6.645E-03 
9.911E-03 
1.951E-03 
3.376E-03 
5.843E-03 
9.644E-03 
1.194E-02 
1.319E-02 
1.359E-02 
1.349E-02 
1.320E-02 
1.280E-02 
1.207E-02 



Table C6. Parameters for SMC (D) 



A (urn) 


a 


P 


7 


0.100 


1.043E+00 


-2.622E-01 


1.688E-02 


0.114 


6.491E-01 


-1.392E-01 


5.393E-03 


0.129 


3.831E-01 


-6.501E-02 


-1.364E-03 


0.143 


2.148E-01 


-2.229E-02 


-4.045E-03 


0.157 


-5.989E-02 


3.263E-02 


-4.722E-03 


0.171 


-3.110E-01 


5.495E-02 


5.624E-04 


0.186 


-4.910E-01 


4.419E-02 


9.322E-03 


0.200 


-5.982E-01 


2.121E-02 


1.698E-02 


0.214 


-6.509E-01 


5.248E-03 


2.060E-02 


0.229 


-6.788E-01 


-9.267E-03 


2.328E-02 


0.243 


-6.912E-01 


-3.038E-02 


2.766E-02 


0.257 


-7.093E-01 


-3.737E-02 


2.829E-02 


0.271 


-7.358E-01 


-2.737E-02 


2.429E-02 


0.286 


-7.384E-01 


-3.690E-02 


2.616E-02 


0.300 


-7.807E-01 


-4.948E-03 


1.573E-02 


0.378 


-8.308E-01 


-1.032E-03 


1.058E-02 


0.475 


-8.621E-01 


-4.739E-03 


8.648E-03 


0.599 


-8.740E-01 


-1.698E-02 


1.029E-02 


0.754 


-8.810E-01 


-2.857E-02 


1.204E-02 


0.949 


-8.860E-01 


-3.811E-02 


1.341E-02 


1.194 


-8.912E-01 


-4.398E-02 


1.398E-02 


1.504 


-8.970E-01 


-4.654E-02 


1.390E-02 


1.893 


-9.034E-01 


-4.678E-02 


1.339E-02 


2.383 


-9.100E-01 


-4.562E-02 


1.263E-02 


3.000 


-9.166E-01 


-4.365E-02 


1.175E-02 



